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ABSTRACT 

Using 766 compact objects found in a systematic survey of the galactic Plane in the 1612-MHz 
masing OH line, new light is cast on the infrared properties of evolved stars on the thermally-pulsing 
asymptotic giant branch and beyond. The usual mid-infrared selection criteria for post-AGB, based 
on IRAS colours, largely fail to distinguish early post-AGB stages. A two-colour diagram from much 
narrower-band MSX flux densities, with bimodal distributions, provides a better tool to do the latter. 
Four mutually consistent selection criteria for OH-masing red proto-planetary nebulae are given, as well 
as two for early post-AGB masers and one for all post-AGB masers including the earliest ones. All 
these criteria miss a group of blue, high-outflow post-AGB sources with 60-/im excess; these will be 
discussed in detail in Paper II. The majority of post-AGB sources show regular double-peaked spectra 
in the OH 1612-MHz line, with fairly low outflow velocities, although the fractions of single peaks and 
irregular spectra may vary with age and mass. The OH flux density shows a fairly regular relation with 
the stellar flux and the envelope optical depth, with the maser efficiency increasing with IRAS colour 
i?2i- The OH flux density is linearly correlated with the 60-/mi flux density. 
Subject headings: Stars:AGB and post-AGB 



1. INTRODUCTION 

Maser emission in the satellite line of groundstate OH 
at 1612.231 MHz occurs in a variety of stellar objects. 
Star-forming regions (SFR), super giants, visible (Miras) 
and invisible (OH/IR) oxygen-rich stars at the tip of the 
asymptotic giant branch (AGB) and post-AGB transi- 
tion objects or even young planetary nebulae can all show 
maser emission at 1612 MHz. Spectral shapes and variabil- 
ity are thought to differ for the different types of object 
(see eg. te Lintel Hekkert & Chapman 1996). 

Objects are often selected as AGB or post-AGB candi- 
dates by the two colours from their flux densities at 12 /im, 
25 fim and 60 /im measured by IRAS (InfraRed Astronom- 
ical Satellite), (eg. te Lintel Hekkert et al. 1991; Lewis 1994 
and references therein). In the resulting two-colour dia- 
gram, a "main sequence" for oxygen-rich AGB stars can 
be identified, where the location of an object is determined 
by its age, initial mass and possibly metallicity (see eg. Be- 
dijn 1987; Likkel 1989). This sequence is called the "evo- 
lutionary" sequence (van der Veen & Habing 1988). To 
test the usual selections, especially of post-AGB objects, 
and potentially find less ambiguous criteria, in this paper 
the approach is, for the first time, from an OH-selected 
sample, finding its IR properties, rather than finding the 
OH properties of an IR-selected sample. Moreover, IRAS 
colours are combined with colours at shorter infrared wave- 
lengths from the Midcourse Space Experiment (MSX) and 
2-Micron All Sky Survey (2MASS). 

The sample consists of 766 compact, OH-masing 
sources 2 in the galactic Plane (for all details see Sev- 
enster et al. 1997a,1997b,2001). The IRAS (12,25,60,100 
/im), MSX (4,8,12,15,21 fim) and 2MASS identifications 
(J,H,K) were taken from the respective public databases. 



Matches were found for 587 (IRAS, within 3c error el- 
lipse), 494 (MSX, within 5" ) and 194 (2MASS, within 
3" ) objects, respectively. Note that the 2MASS mission 
has not yet covered the entire section of the plane cov- 
ered by the OH sample. For many sources that do have a 
valid positional identification, several or most of the flux 
densities may not have been measured properly. For in- 
stance, out of the 587 IRAS identifications, only 240 have 
well-determined 12-/im , 25-/im and 60-/im flux densities 
(quality flag 2 or 3). The relatively small percentage of 
IRAS identifications is caused largely by the confusion at 
the very low latitudes of our survey (\b\ < 3°. 25). For the 
same reason, the flux-density measurements, especially at 
60 fim. and 100 fim, are likely to be somewhat overesti- 
mated even if they are not flagged as "upper limits" in the 
IRAS data base. With the exception of a few likely su- 
pergiants, none of the sources could be identified with an 
optical counterpart. The combination of low latitude and 
OH-based selection means that lowest-mass AGB stars 
(< 1 Mq) may be underrepresented in the sample. 

The aim of this analysis is to find generic differences 
between post-AGB and AGB objects, both in infrared 
and maser properties, and signs of post-AGB evolution. 
We will define a post-AGB stage as the phase after the 
last thermal pulse, when variability has ceased. The term 
proto-planetary nebula (PPN) will indicate a slightly later 
stage, when the mass-loss rate has dropped by orders of 
magnitude and the near-infrared or even optical emis- 
sion is becoming much stronger with respect to the mid- 
infrared. Sometimes such objects are called OHPN, as 
they may already show radio continuum (PN), but still 
harbour masers (OH). However, as long as an object shows 
OH-maser emission, especially at 1612 MHz as our sample 
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Fig. 1. — The histograms of the two main IRAS colours (R21 = [12 
well-determined flux-density measurements for all three bands involved. 




R32 



25] and R32 = [25 — 60]) for all OH-selccted sources with 



does, it is per definition not a full-blown PN (Zijlstra et 
al. 1989; Lewis 1989). 

In Section 2, we will introduce an MSX two-colour di- 
agram. In Section 3, we discuss the role of the OH out- 
flow velocity and the spectral shapes of different types of 
objects and in Section 4 the near-infrared properties are 
added to the discussion. In Section 5 some selection cri- 
teria for late versus early post-AGB objects are summa- 
rized and evidence from known sources presented, then in 
Section 6 a new expression for OH 1612-MHz maser effi- 
ciency in (post-) AGB objects is given. Conclusions are 
presented in Section 7. 

2. DISTRIBUTIONS OF MID-INFRARED COLOURS 

All mid infrared (MIR) colours are defined as [a — b] = 
2.51og(Sb/Sa) W1 th S flux density in Jy and a,b wave- 
length in /xm. For IRAS, the usual names R21 = [12 — 25] 
and i?32 = [25 — 60] are used, as well as R43 = [60 — 100]. 
The 12-/xm bands of the two satellites are different : for 
MSX it covers 11.1 /j,m to 13.2 /j,m, for IRAS 8.5 /um to 
15.5 /im. This means, most importantly, that the 9.7-/zm 
silicate feature is in the MSX 8-/^m band (6.8 /xm to 10.8 
/im) whereas it is in the IRAS 12-^tm band. Only in very 
cold material may the silicate feature shift into the MSX 
12-/zm band (Zhang & Kwok 1990). The MSX 15-^m 
band is almost entirely included in the IRAS 12-/zm band 
and similarly the MSX 21-/im band is almost entirely in- 
cluded in the IRAS 25-/im band. 

In Fig. 1&2, histograms are shown for five MIR colours. 
The distributions for [8 — 12] and [15 — 21] immediately 
stand out as curiously bimodal when compared to the 
other three, with divisions at ~0.9 and ~0.7, respectively. 
To follow this up, we use those two colours for the MSX 
two-colour diagram. All sources with these MSX colours 
well-defined are plotted in Fig. 3. The bulk of the sources 
is in the bluest quadrant, with little spread in colour, but 
a second group is at much redder [8 — 12] with a large 
spread in [15 — 21]. 

In Fig. 4, the traditional IRAS colour-colour diagram 
is shown, with the evolutionary sequence and the regions 
as laid out by van der Veen & Habing (1988) and sym- 
bols defined by MSX colour. The circles are found in the 
region where one typically expects late post-AGB objects 
by IRAS colours (region V and further; Bedijn 1987; van 
der Veen & Habing 1988 ("non-variable OH/IR stars"))- 



The bulk of the main MSX group is along the evolution- 
ary sequence, except for a surprisingly large group in the 
left part of region VIII and above. The crosses are found 
mostly straddling the border between V and VIII, below 
the "bipolar outflow" region, defined recently by Zijlstra et 
al. (2001). A solitary '+' is right on the edge of where one 
expects to find OH-emitting SFRs (Braz & Sivagnanam 
1987); in this region there is strong overlap between PPN 
and SFRs. 

The bimodality of the two MSX-colour distributions 
suggests the existence of at least two intrinsically different 
types of objects in the sample. Combining this with what 
one already understands from the traditional IRAS colour- 
colour diagram (van der Veen & Habing 1988), it seems 
that the transition from blue (<0.9) to red (>0.9) [8 - 12] 
coincides with the transition off the AGB. The transition 
from blue (<0.7) to red (>0.7) [15 - 21] could well indi- 
cate a next evolutionary transition, further away from the 
evolutionary sequence (Fig. 4), such as the drop in mass- 
loss rate and start of the fast PPN wind (cf. van Hoof et 
al. 1997). Hence, the MSX diagram possibly distinguishes 
between different post-AGB stages and provides an eas- 
ier tool to separate post-AGB stars from AGB stars, but 
does not distinguish between different AGB objects like 
the IRAS diagram with its evolutionary sequence. The 
two diagrams should be used in combination. 

3. THE ROLE OF THE OH-SHELL OUTFLOW VELOCITY 

The outflow velocity V cxp of a typical OH maser spec- 
trum is thought to be related to the luminosity (age, initial 
mass) and the metallicity (age, galactic coordinates) of the 
central star, via 

L. = V c %Z-z (1) 

according to van der Veen (1989; see also Habing et 
al. 1994). Not all OH maser spectra are regular and 
double-peaked, so the outflow velocity is not always well- 
defined. Specifically, some suspected post-AGB objects 
may have spectra ranging over almost 100 kins -1 (e.g. 
spectrum 2 in Fig. A7; see also te Lintel Hckkert & Chap- 
man 1996), but be listed as single-peaked with just the 
peak velocity, whereas similar objects do have an asso- 
ciated outflow velocity (Fig. Al(4 th -last)). Such single- 
peaked objects are intrinsically different from star-forming 
regions, that commonly - especially with sufficient spatial 
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Fig. 2. — The histograms for the three main MSX colours, for all sources with well-determined flux densities in both bands of each colour. 
The [8 — 12] and [15 — 21] distributions look very different from the [12 — 15] distribution, as well as from the distributions in Fig. 1, and seem 
bimodal with minima around 0.9 and 0.7, respectively. 
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Fig. 3. — A two— colour diagram for MSX colours, split into four quadrants according to the apparent minima in the bimodal distributions 
shown in Fig. 2. Objects are plotted with different symbols accordingly (circles in the first quadrant, plusses in the second, etc.), that will be 
used to relate their locations in other diagrams to their MSX colours. The dashed line shows the colours for black-body radiation at various 
temperatures (tags). 
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Fig. 4. — The two-colour diagram for IRAS colours, with symbols as defined by the quadrants in Fig. 3. The solid curve indicates the 
evolutionary sequence and the dashed boxes the regions as defined by van der Veen & Habing (1988). The dashed line separates roughly the 
red corner of OH-emitting SFRs (Braz & Sivagnanam 1987; see Section 5.3 ). (Non-OH-emitting SFRs would be found to the left of this 
region at R32 >2.) The dotted rectangle is the region classified as "bipolar outflow" region by Zijlstra et al. (2001). Note that compared to 
Fig. 3 there are far fewer objects, as now the IRAS colours, as well as the MSX colours, have to be well-defined as well as the MSX colours. 
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Fig. 5. — The MSX [8 — 12] colour plotted against OH outflow velocity. The general group of AGB sources with [8 — 12] ~0.4 ranges over 
all outflow velocities. There is a discrete group of sources with much redder colour, indicating they are post-AGB and PPN sources. Their 
outflow velocities arc limited to the range 9-15 kms" 1 . The same trend is seen - less clearly - for R21 and [15 — 21], but not for any of the 
other IR colours. 




Fig. 6. — As Fig. 4, with symbols according to outflow velocity. The extremely red single-peaked sources ('+') in the top panel are 
most likely to be SFRs. The PPN region is populated mostly by sources with outflow velocities below 15 kms -1 (as derives from Fig. 4&5). 
Interesting is the region to the left of the evolutionary sequence, with < R32 < 2, where sources have mostly outflow velocities higher than 
15 kms" 1 . Before (Fig. 3& 4), these sources were found to have very blue MSX colours. Irregular sources (Section 3) are found in and around 
the bipolar-outflow region (lower panel). The spectra of the irregular sources are shown in Fig. A7. 
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resolution - have just one very narrow bright peak in their 
1612-MHz spectrum, single-peaked PN (Fig. Al(13)) or 
single-detected AGB stars, that simply have one of the 
peaks missing due to detection limits (cf. Fig. Al(15)). 
Hence, selecting K xp =0kms _1 results in a hotchpotch of 
objects. 

We want to find out what OH outflow velocity and spec- 
tral "shape" tell us about the evolutionary stage of an ob- 
ject. In Fig. 5, [8—12] is plotted against outflow velocity. A 
population of red sources is clearly separated from the rest, 
with outflow velocities between 9kms~ 1 and 15kms _1 . 
The separation is at [8 - 12] -0.9, as before (Fig. 2&3); 
for [15 — 21] a similar but less clear separation is found at 
—0.7 . Thus, interestingly, in the first and fourth quadrant 
of Fig. 3, all double-peaked sources are found to be in this 
very limited range of outflow velocities. 

In the IRAS colour-colour diagram (Fig. 6a), region V 
and beyond is populated mainly by lower-outflow veloc- 
ity objects, as follows from the previous two figures. The 
single-peaked sources in the reddest corner of this diagram 
are most likely SFRs. On the evolutionary sequence, the 
spread in outflow velocities is largest. We do not find 
the increase of outflow velocity with i?2i for blue sources 
that is reported in Chen et al. (2001), but our OH-selected 
sample is probably too red to see this. The higher outflow 
velocities prevail in the left-hand sides of regions IV and 
VIII, where we noted relatively high numbers of sources 
before (Section 2). These sources may also be post-AGB 
sources (van der Veen & Habing 1988); they will be dis- 
cussed in Paper II (Sevenster 2002). 

Post-AGB stars, typically selected from region V, are 
often connected to irregular, extreme-velocity outflows (te 
Lintel Hckkert & Chapman 1996). However, looking at the 
OH spectra for these sources (Fig. Al), they correspond to 
OH-selected sources with mostly normal, double-peaked, 
intermediate -outflow spectra (Table 1). The same is true 
for the higher - outflow sources to the left of the evolution- 
ary sequence, although the fraction of irregular spectra is 
higher (Fig. A5, Table 1). 

This is not to say that there is no subset of irregular 
sources with reasonably confined properties. In Fig. A7, 
the spectra are shown for all sources from the full sam- 
ple that have velocity widths (W20) much larger than the 
peak separation (AV = 2V cxp ). This criterion is used 
as an indicator of non-spherical, irregular outflow. Their 
IRAS colours (Fig. 6b) clearly indicate a connection to the 
bipolar-outflow region of Zijlstra ct al. (2001). The actual 
velocity widths are up to 70 kins" 1 . 

In summary, although very irregular OH spectra may 
always indicate a post-AGB host, the reverse may not be 
true. Unless the strong reddening in [8 — 12] and i?2i is 
unrelated to the onset of post-AGB processes, OH 1612- 
MHz spectra of post-AGB objects seem to be mainly reg- 
ular and show no indication of irregular mass loss or in- 
creased outflow velocity, spherical nor bipolar. For the 
bulk of the red post-AGB objects, the outflow velocities 
are even lower than average: between 9 kms -1 and 15 
kms" 1 . Most likely, this indicates that the OH masers in 
general disappear before the changes of the central wind 
truly influence the envelope at such large radii. 

4. NEAR-INFRARED PROPERTIES 



From the 2MASS database, J (1.25 /an), H (1.65 /an) 
and/or K (2.17 /an) magnitudes were obtained for 194 
sources. Colours are used as uncorrected J-H and H-K. 
The two-colour diagram is shown in Fig. 7. Using the 
symbols from Fig. 3, post-AGB stars are where one may 
expect them (e.g. Garcia-Lario et al. 1997) but the dia- 
gram can hardly be used as a diagnostic for any particular 
type of object (see also Fig. 10c). 

An alternative to the IRAS colour-colour diagram was 
proposed by van Hoof et al. (1997), using [12] — [25] 
and K-[12] (Fig. 8), with IRAS magnitudes [12] = 
2.51og(59.5/5i 2 ) and [25] = 2.51og(13.4/5 25 ) following 
Oudmaijer et al. (1992). The PPN symbols from Fig. 3 
(circles) are exclusively at [12] — [25]<;3.5. In the bottom 
panel the selection is based on OH properties, as well as 
[15 — 21] . The latter colour is used for the selection, rather 
than [8 — 12], because it may single out the further-evolved 
post-AGB sources (i.e. PPNe). Selecting 14 xp between 
9kms~ 1 and 15kms _1 (Fig. 5) and the sum of the OH 
flux densities of the two peaks larger than 1 Jy, we find 
a sample (circles) that occupies the same region as the 
circles in Fig. 8(a). The dashed line in these diagrams is 
drawn to separate the circles from the other objects. 

The post-AGB selection of Fig. 8(b), based largely on 
OH spectral properties, gives without further assumptions 
[12] - [25] > 3 and K-[12] < 9, as well as [8 - 12] > 0.9. 
This once more confirms that the sources in the first quad- 
rant of Fig. 3 are "standard" post-AGB sources, moving 
to higher R21 in the IRAS two-colour diagram (Bcdijn 
1987; van der Veen & Habing 1988), and that these are 
mostly regular and double-peaked (Fig. A3). A second 
mixed near-mid-infrared two-colour diagram is shown in 
Fig. 9. The red post-AGB group separates clearly from 
the main oxygen-rich colour sequence (see van Loon ct 
al. (1998)). However, some objects that are still in the 
third quadrant of the MSX diagram - and on the IRAS 
evolutionary sequence - are also above the dashed curve 
and quite distinctly separated from the main colour se- 
quence. This may indicate that they are in the very first 
stages of post-AGB colour change. 

In conclusion, the NIR magnitudes of objects are useful 
to distinguish (post-) AGB evolutionary stages primar- 
ily in combination with the 12-/an and 25-/an IRAS flux 
densities. 

5. OH SPECTRA OF IR SELECTED SOURCES 

5.1. Selecting post-AGB sources 

Summarizing the previous three sections, standard 
"cold" evolved OH PPN objects can thus be selected vir- 
tually consistently in four ways. 

• (1) The traditional IRAS selection : i? 32 ^1.5 and 
i?2i^l-4 (Fig. 4). Spectra are shown in Fig. Al. 

• (2) From MSX colours : [8-12] > 0.9 and [15-21] > 
0.7 (Fig. 3). Spectra are shown in Fig. A2. 

• (3) From OH 1612-MHz properties : V cxp = 12 ± 
3kms _1 and S'oh blue + Soh red > 1 Jy (resolution 
^kms- 1 ), with [15 - 21] > 6.7 (Fig. 8b). Spectra 
are shown in Fig. A3. 

• (4) With NIR : [12] - [25] > 4.164 - 0.164 (K - [12]) 
and K— [12] < 9 (Fig. 8a). This selection is mostly a 
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Fig. 7. — The 2MASS two-colour diagram (Section 4) with symbols defined by MSX colours (Fig. 3). The dashed line connects black-body 
colours for several temperatures. For J— H < 0.4, Garcia-Lario et al. (1997) find exclusively planetary nebulae (PNe can be found at redder 
colours). None of our OH-masing sources are that blue, which would agree with the idea that masers, especially the 1612-MHz OH maser, 
cannot exist in the environment of a truely ionized, high-wind nebula. 
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Fig. 8. — The two-colour diagram using combined MIR and NIR colours as defined in Section 4, following van Hoof ct al. (1997). a. 
Symbols defined by MSX colour (Fig. 3). b. Symbols defined mainly by OH properties. The circles indicate sources with outflow velocities 
between 9kms _1 and 15kms _1 , the sum of the two peak-flux densities larger than 1 Jy and [15 — 21]>0.7. They also have [15 — 21] >0.7. 
Based on our previous findings, these sources are likely "standard" PPN sources. Note that all circles have K— [12] < 9, which is relevant in 
the context of Fig. 10. 




Fig. 9. — The two— colour diagram using combined MIR and NIR colours as defined in Section 4, following van Loon et al. (1998). The solid 
curve indicates their oxygen-rich AGB sequence; the dashed curve is our separation between AGB and post-AGB objects. The post— AGB 
sources, from quadrants I and IV (Fig. 3), are all found above this line. Quite a few sources from quadrant III are found in the post-AGB 
region, however, which suggest that this diagram might provide a tool to distinguish the earliest post-AGB stars. Star-forming regions are 
found at K-[12] >9 (Fig. 10). 
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R21 K-[12] H-K H-K 

Fig. 10. — Two— colour diagrams for a sample of known OH— masing SFRs (Testi et al. 1998). Circles indicate reliable flux-density mea- 
surements and crosses colours from flux-density upper limits. In a, the IRAS two— colour diagram is plotted, showing the "OH-masing SFR" 
region (Section 5.3). In b, the near-mid— infrared diagram shows that SFRs all have K— [12] > 9 but are above the post— AGB line; compare to 
Fig. 8. Hence, the region 9<K-[12]<12 is shared by post-AGB stars and SFRs, but SFRs do not venture into the "PPN" region (Fig. 8b). 
In c, SFRs are seen to are spread over all NIR colours, hence a good fraction has H— K higher than post-AGB objects. In d, this means SFRs 
and post-AGB sources overlap only in the region K— [12] > 9 and H— K < 1.5. 
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Fig. 11. 
with X21 : 
ratio. 



- According to van der Veen (1989), the OH luminosity for objects along the evolutionary sequence is proportional to L\ X^, 
525/512. From this log-log plot, slightly different powers are found, 1.5 for the stellar luminosity and 3.5 for the flux-density 



doubles, al 

o : R21 > 1.5 
x : R21 > 0.5 
« : R21 < 0.5 
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Fig. 12. — Same as Fig. 11, including (double-peaked) sources away from the evolutionary sequence. The power for X21 is 2 in this case. 
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Fig. 13. — The maser "efficiency", SoH/Sboi (Section 6), is plotted logarithmically against IRAS colour R21. For the bluest sources, it 
is <x X^, for the reddest sources oc (both Son and S bol are almost constant). Again, only double-peaked sources arc plotted to avoid 
SFRs. 



combination of the previous ones, but it can be used 
for sources with known K magnitude and no MSX 
or 60-/um association. 

As seen in all the corresponding figures, the selection cri- 
teria, partly overlapping, are all entirely compatible, with 
just one or two sources breaking the rules. 

Early post-AGB sources are selected by : 

• (5) [15 - 21] < 0.7 and [8 - 12] > 0.9 (Fig. 3,4&8(a)), 
spectra in Fig. A4 

• (6) [12] - [25] > 4.164 - 0.164 (K - [12]) and 
K— [12] >9 (Fig. 8a); these colours are shared by 
SFRs (Fig. 10) 



Post-AGB sources in general, including the very earliest, 
may be selected by : 

• (7) (H-K) + 0.05 < 0.03 ((K-[12])+0.807) 2 , where 
early and late post-AGB objects can be distin- 
guished by [15 — 21] smaller or larger than 0.7 or 
even by [8 — 12] smaller or larger than 0.9, respec- 
tively (Fig. 9, Fig. 3). 



A group of blue possible post-AGB stars is found at : 

• (8) i? 32 > -0.2 and R 2 i > 0.2, to the left of the evolu- 
tionary sequence (i?32 > —2.15+0.35 exp(1.5 i?2i))- 

These are probably high-mass post-AGB sources (see Pa- 
per II)). Non-spherical-outflow objects (Fig. A7&6b) are 
likely to be found roughly in the region defined by Zijlstra 
et al. (2001). 

5.2. Known PNe in the sample 

Taking all the sources from our OH sample that are as- 
sociated with a "Planetary Nebula" within 5" according 
to SIMBAD, 10 out of 16 do satisfy the OH criteria of 
outflow velocity and flux density. Four of the 16 have a 
2MASS association, with K— [12] < 9.4 and above the lines 
in Fig. 8&9. 



Fourteen have an IRAS association, with 8 in the tra- 
ditional PPN region to the right of the evolutionary se- 
quence and 1 to the left, with V exp > 15 kms" 1 and very 
blue MSX colours. Of the 12, out of 16, sources with MSX 
association, 5 are in the first quadrant (Fig. 3). 

Spectra of all 16 sources are given in Fig. A6. Even 
of this group of selected "Planetary Nebulae", 75% 
have regular, double-peaked spectra. Even two known 
"bipolar planetary nebulae" are in this group of 16 : 
OH353.844+02.984 (Fig. A6(3); Kwok et al. 1996) and 
OH348.813-02.840 (Fig. A6(6); Hrivnak et al. 1999). The 
latter satisfies all our selection criteria for which its prop- 
erties are known, including slow, regular outflow. 

In summary, sources classified as "planetary nebulae" in 
fact cover a range of post-AGB evolutionary stages. As 
we noted in the caption of Fig. 7, a source that still har- 
bours 1612-MHz emission is unlikely to be a full-blown 
PN. It would be good practice to name such objects PPN 
or OHPN (see Introduction). 



5.3. Colours of SFR sources 

It is hard to define clear-cut selection criteria to find 
SFRs in an OH sample, especially to tell them apart from 
PPN. A sample of known OH-masing SFRs with NIR 
data, from the literature (Testi et al. 1998), is used to check 
their infrared behaviour with respect to post-AGB stars. 
In Fig. 10, several two-colour diagrams are given for com- 
parison with the diagrams in the rest of the paper. The 
dotted line in Fig. 10(a) outlines the OH-masing SFR re- 
gion, slightly adapted from Braz & Sivagnanam (1987), 
who claim that non-OH-masing SFRs are found to the 
left of this region, at lower i? 2 i but similar i? 32 . 

Interestingly, in K-[12] versus [12] - [25] (Fig. 10b), the 
SFRs lie above the empirical line drawn in Fig. 8, but at 
higher K— [12] than the PPN sources. The selected sources 
in Fig. 8b and the SFRs in Fig. 10b are perfectly separated, 
although there are some earlier post-AGB sources in the 
region K-[12]>9 (Fig. 8a). In Fig. lOd, the SFRs are 
also in the post-AGB region, again at higher K— [12] and 
mostly higher H— K (Fig. 10c). Three SFRS with a reli- 
able MSX identification with (some) well-determined flux 
densities are all in quadrant II. 
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Table 1 

Fractions (in %) of double (D), single (S) and irregular (I) spectra for several post-AGB selections 
plus comparison PN and AGB samples (judged by eye from the figures in appendix A). The first four 

ROWS CORRESPOND TO SELECTIONS 1,8,2 AND 5 IN SECTION 5.1, RESPECTIVELY 
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comment 


Right 
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15 


iras (low-mass) 


Left 


65 
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iras (high-mass) 
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22 


msx (old) 


Q4 


75 
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msx (young) 


PN 
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10 


15 


simbad (Section 5.2) 


AGB 


80 
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15 


iras (on evol. seq.) 



* * **■■ S *% t 
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R21 



Fig. 14. — The ratio of OH-pcak flux density and 60-fim flux density plotted logarithmically against IRAS colour R21 for all double-peaked 
sources. Of all IRAS flux densities, S60 correlates most strongly, almost linearly, with Soh- The 560 contains the 53-fim pumping line. The 
dashed line is the relation found by Chen ct al. (2001) for S011/S35, with S35 the interpolated flux density at the 35-^tm pumping line. 
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6. MASER-INFRARED CORRELATION 

In order to determine a maser "efficiency" in terms of 
OH flux density as a function of stellar flux, we derive the 
bolometric flux, Sboi = £12 * BC12, using the bolometric 
correction 

BC 12 = 0.7 + 2.9 cxp(-3.0 R 21 ) + 0.9 cxp(0.7 R 21 ) (2) 

from van der Veen & Breukers (1989). If all the stellar 
light is recycled in the infrared, then S* = S\> a \- Van 
der Veen (1989) gives, for sources on the evolutionary se- 
quence, Soh = S 1 * X$i, where X 21 = S 2 ^/S\ 2 . We find 
that this relation is not exactly right and also that the 
scatter around it is very large. Rather, for sources close 
to the evolutionary sequence (Fig. 11), 

Sob. = Sbai X 2i (3a) 
and for the whole sample in general (Fig. 12) 

Sob. = Sboi-^ii- (3k) 

However, the "maser efficiency", defined as Soh/SW, is 
proportional to for the very bluest sources (Fig. 13), 
since of course S'boi itself varies strongly with X 2 i (Equa- 
tion 2). For the reddest sources, the maser efficiency is 
virtually constant. With a scatter of about a decade in 
these log-log plots these relations cannot be used on a 
one-by-one basis. We can conclude, however, that sta- 
tistically the OH flux density of circumstellar masers is a 
well-behaved function of "pumping" (Sboi) an d "optical 
depth" (i? 21 ). 

Moreover, there is a very clear, almost linear correlation 
between the 60-/zm flux density and the OH flux density 
(Fig. 14). The 60-^m band contains the 53-/xm pump- 
ing line, which may explain the correlation with "available 
pumping photons" , even though the line would be too nar- 
row to influence the wide-band flux. The red post-AGB 
sources have surprisingly similar flux-density ratio to the 
evolutionary-sequence sources. In Fig. 14, we plot the 
slope of the relation found by Chen et al. (2001) for the 
ratio S0H/S35, with the IRAS flux density at the main 



pumping line at 35 /im derived by interpolation. Not only 
is their relation not linear, it is also less tight, so it seems 
the Seo is a more uniform indicator of pumping than the 
interpolated S35, possibly due to contaminating lines in 
the S25 band. 

7. CONCLUSIONS 

An extensive comparison is given between several in- 
frared and OH-maser properties of an OH-selected sample 
of objects in the galactic plane. All the OH specifics in this 
paper apply only to the 1612-MHz maser line and may be 
(quite) different when studying the main lines. Four differ- 
ent selection criteria are presented for finding red proto- 
planetary nebulae with OH masers, that all give consistent 
selections. Most efficient in separating post-AGB stars is 
the MSX two-colour diagram for [8 - 12] and [15 - 21] 
that splits into four clear quadrants, containing primar- 
ily PPNe, SFRs, AGB objects and post-AGB stars, re- 
spectively. Two selection criteria are presented for finding 
early post-AGB OH stars and one for post-AGB stars of 
all ages. A significant group of blue (likely) post-AGB 
stars with 60-/xm excess remains unselected by all of the 
above criteria; these will be discussed in Paper II. 

None of the selected post-AGB groups show any 
true prevalence of irregular or extreme-velocity outflows, 
rather about 70% of their OH spectra are typical of regu- 
lar, spherical, thin-shell outflow with velocities between 9 
kms" 1 and 15 kms -1 . However, the sources with very ir- 
regular spectra are located roughly in the bipolar-outflow 
IRAS region defined by Zijlstra et al. (2001). The reddest 
PPN may have a slightly higher fraction of genuine single- 
peaked sources. 

Maser efficiency, defined as the ratio of OH luminos- 
ity to bolometric luminosity, increases steeply with IRAS 
colour i?2i for blue sources, but is almost constant for the 
reddest sources. Of all IRAS flux densities, the OH maser 
correlates most strongly, and linearly, with the 60-fim flux 
density. 

MS thanks the Leidse Sterrewacht for kindly providing 
facilities to finish the final version of this paper. 
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APPENDIX 
SPECTRA 






A^AA|f 




Fig. Al. — All sources with R32 <1.5 and R21 >1.4 (item 1 in Section 5.1). Note that all spectra in this appendix have flux densities that 
are NOT corrected for primary-beam dilution. In the text, quoted flux densities ARE corrected. 
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Fig. A3. — All sources with [15-21]>0.7 and 9kms V <15 kms and sum of the blue— and red— shifted OH peak flux densities larger 
than 1 Jy (item 3 in Section 5.1, Fig. 8). This gives [12]-[25]>3 and [8-12]>0.9 . 
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Fig. A4.— Sources with [15-21]<0.7 and [8-12]>0.9 (item 5 in Section 5.1). Second last source has [15-21]~0.7 though 
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Fig. A5. — All sources with R32 > —0.2 and .R21 >0.2 to the left of the evolutionary sequence in the IRAS two-colour diagram (item 8 in 
Section 5.1, Fig. 4). 
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Fig. A8. — Thermally-pulsing AGB sources, on the IRAS evolutionary sequence (solid curve in Fig. 4). 



